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We present results from numerical simulations of full, two flavor QCD thermodynamics at A^i = 4 with domain 
wall fermions. For the first time a numerical simulation of the full QCD phase transition displays a low temperature 
phase with spontaneous chiral symmetry breaking but intact flavor symmetry and a high temperature phase with 
the fuU 517(2) X SU{2) chiral flavor symmetry. 



1. Introduction 

QCD thermodynamics has been extensively 
studied using numerical simulations with stag- 
gered or Wilson fermions. However, both for- 
mulations break the chiral symmetry of the the- 
ory. The symmetry is recovered together with the 
Lorentz symmetry as the continuum limit is ap- 
proached. In the past few years a novel fermion 
regulator was developed that provides a way of 
controlling the amount of chiral symmetry break- 
ing at any lattice spacing. Domain wall fermions 
utilize an extra space-time dimension with free 
boundary conditions to separate the two chiral 
components of the Dirac spinor. The components 
are localized on the opposite boundaries (walls). 
If the extent of the fifth "dimension" , Lg, is infi- 
nite the two chiral components decouple and the 
theory has the full chiral symmetry, even at finite 
lattice spacing. For practical numerical simula- 
tions, Lg must be finite and as a result there is a 
small mixing of the chiral components resulting in 
a residual mass. The important point is that the 
size of this residual chiral symmetry breaking can 
be controlled at any lattice spacing by increasing 
Lg. For the first time the approach to the chiral 
limit has been separated from the approach to the 
continuum limit. 

In free field theory the localization is exponen- 
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tial and the effective quark mass is given by 



(free) 
^eff 



mo(2 - mo) [ruf + (1 - mo) 



(1) 



where < mo < 2 is the domain wall "height" 
and m^ an explicit fermion mass. In the interact- 
ing theory there is numerical non-perturbative ev- 
idence that the general features of eq. | still hold 
with mo being renormalized ||], |||. 

Another unique property of domain wall 
fermions is that the Lg = oo limit can be stud- 
ied using the overlap formalism In that 
limit it can be shown that for topological gauge 
field backgrounds there are exact fermionic zero 
modes. It has been demonstrated that for masses 
in the range of interest the effects of finite Lg on 
the zero modes can be made arbitrarily small for 
classical Q, ||] and for quantum [Q, |^ gauge 
field backgrounds. Domain wall fermions possess 
the key ingredient in reproducing anomalous ef- 
fects. 

For more details on domain wall fermions see 
the review [|l0| and references within. 

2. Numerical Simulations 

Dynamical domain wall simulations have been 
done in the past for the Schwinger model in 
the Lg — oo limit (overlap formalism) |pT[| and 
for finite Lg [|j using standard hybrid Monte 
Carlo (HMC) techniques. In this work we 
simulated QCD using standard HMC at finite 
Lg. The gauge fields are defined on the four- 
dimensional lattice while the fermion fields on 
the five-dimensional one. The five-dimensional 
Dirac operator is as in |l3] with even-odd pre- 
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conditioning. The bulk effects of the Lg heavy 
flavors were subtracted [H by introducing five- 
dimensional bosonic fields as in §. We used the 
<I> algorithm with trajectory length of 0.5 and step 
sizes ^ 0.01 — 0.02 resulting in acceptance rates 
^ 80%. The standard conjugate gradient inver- 
sion algorithm (CG) was used with a typical num- 
ber of CG iterations ranging from 50 — 200. We 
used initial configurations in the opposite phase 
and 200 — 800 thermalization sweeps. The com- 
putational cost is linear in Lg. 

3. The Phase Transition 

In order to investigate the feasibility of study- 
ing QCD thermodynamics with domain wall 
fermions, to locate the critical coupling, and to 
investigate the parameter space of this new regu- 
lator {mo,mf, Lg) we performed a large number 
of simulations at various couplings on 8^ x 4 and 
IG'^ X 4 lattices using the QCDSP machine. This 
was possible due to the robustness of the machine 
and software that allowed us to split a 200 Gflops 
portion of it to 7 independent, 6 Gflops machines 
and 6 independent, 25 Gflops machines. 

A sharp change on the value of the chiral con- 
densate (tjj'ip) and the magnitude of the Wilson 
line I IV I was observed as (3 = 6/gQ (with go being 
the gauge coupling constant) was varied. This 
can be seen in figure 1 for an 8"^ x 4 lattice with 
Too = 1.9, rrif = 0.1 and = 12. The fits 
are to co[ci -I- tanh(c2[/3c — P])]- In the broken 
phase we did a dynamical mass extrapolation of 
{ip'4') at /3 = 5.2 with Ls — 16 and found a non 
zero intercept, 9.1(1) x 10^"^, indicating the ex- 
pected spontaneous breaking of chiral symmetry. 
In the symmetric phase a dynamical mass extrap- 
olation of ('0V') at /? = 5.45 with = 8 re- 
sulted in an intercept of 2.8(1) x 10"'^ but when Lj 
was increased to Lg = 16 the intercept dropped 
to 4.5(1.3) X 10^^ indicating that the full chiral 
symmetry was restored to a high degree. Simi- 
lar results were obtained on a 16* x 4 lattice at 
/3 = 5.45, Too = 1.9 and Ls — 16, resulting in an 
intercept of 4.4(5) x 10~*. The dynamical mass 
data and extrapolations are presented in figure 2. 

The Ls dependence of our results is presented 
in figure 3. The lines are curves of the form 



Co -t- cie~'^^^= that go through the three points. 
The horizontal lines are the (V'V') = cq lines. As 
can be seen in the broken phase (/3 = 5.2) where 
the coupling is larger, the decay rate C2 is larger 
than the one in the symmetric phase (/? = 5.45) 
as expected from the Schwinger model At 
Ls = 16 and /3 = 5.2 ("00) is within 7% of the 
extrapolated value while at /3 = 5.45 it agrees 
within the error bars. 
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Figure 1. 8" x 4, = 12, m/ = 0.1, and 
Too = 1.9. \yV\ (crosses) and (V'V') (diamonds). 
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Figure 2. 8^ x 4 /3 = 5.20 (squares), 8^ x 4 
/3 = 5.45 (diamonds), 16^^ x 4 = 5.45 (crosses), 
with Ls = 16 and toq = 1.9. The stars are the 
TOy = linearly extrapolated values. 

The dependence on toq is more subtle. In free 
field theory toq < corresponds to no light fla- 
vors, < Too < 2 to one and 2 < too < 4 to four 
light flavors. Also, for positive but small Too the 
number of low momentum states becomes small 
(for more details and numerical results see 1^). In 
order to study these effects we located /3c using 
(■0?/') and |iy|, for various values of toq on 8'^ x 4 
lattices. The results are in figure 4. The depen- 
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dence of Pc on tuq can be viewed as an expected 
renormalization effect. 
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Figure 3. 8^ x 4, to/ — 0.1 and too = 1.9. /3 = 
5.20 (crosses), (3 — 5.45 (diamonds). 



Co value suggests that UA{i) remains broken just 
above the transition. 
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Figure 5. 16^ x 4, Ls ^ 16, /3 = 5.45, and 
Too ~ 1.9. The fit is to cq + C2m'j and the star is 
the ruf = extrapolated value. 
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Figure 4. 8"' x 4, = 12 and to/ = 0.1. /3c 
from \W\ (crosses) and from (i/'V') (diamonds). 



5. Conclusions 

We presented results from numerical simula- 
tions of full, two flavor QCD thermodynamics at 
Nt = 4: with domain wall fermions. They indi- 
cate the presence of a low temperature phase with 
spontaneous chiral symmetry breaking but intact 
flavor symmetry and a high temperature phase 
with the full SU{2) x SU{2) chiral flavor symme- 
try. Given this, we interpret the difference seen in 
the TT and 5 screening lengths just above the tran- 
sition as preliminary evidence for the anomalous 
breaking of the Ua{^) symmetry. 



4. The U{1) Axial Symmetry 

Due to the exceptional zero mode properties 
of domain wall fermions the possibly anomalous 
breaking of the Ua{^) symmetry just above the 
transition can be investigated with clarity and 
without subtleties such as zero mode shift effects 
(staggered fermions) or proximity to the Aoki 
phase (Wilson fermions). We measured the dif- 
ference of the screening masses of the 5 and tt as a 
function of to/ on a 16 x 4 lattice at /3 = 5.45 with 
Too = 1.9 and Ls = 16. The data fits to co + C2'iTij 
with x^/dof « 0.5 and cq = 3.1(9) x 10"^^ The 
small /dof and the near-zero, small-mass limit 
of (ipip) indicate that the chiral symmetry break- 
ing effects due to the fermion regulator {Lg < oo) 
are negligible so that this preliminary, non-zero 
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